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Abstract. 

We calculate matter effects on neutrino oscillations relevant for long baseline neutrino oscillation experiments. In partic- 
ular, we compare the results obtained with simplifying approximations for the density profile in the Earth versus results 
obtained with actual density profiles. We study the dependence of the oscillation signals on both ii/Am^,^ and on the 
angles in the leptonic mixing matrix. The results show quantitatively how matter effects can cause significant changes in 
the oscillation signals, relative to vacuum oscillations and can be useful in amplifying these signals and helping one to 
obtain measurements of mixing parameters and the magnitude and sign of Am . 



INTRODUCTION 

In a modern theoretical context, one generally ex- 
pects nonzero neutrino masses and associated lepton mix- 
ing. Experimentally, there has been accumulating evi- 
dence for such masses and mixing. All solar neutrino ex- 
periments (Homestake, Kamiokande, SuperKamiokande, 
SAGE, and GALLEX) show a significant deficit in the 
neutrino fluxes coming from the Sun (1). This deficit can 
be explained by oscillations of the v^'s into other weak 
eigenstate(s), with Am^^, of the order 10"^ eV^ for MSW 
solutions (2) or of the order of 10^'" eV^ for vacuum 
oscillations. Accounting for the data with vacuum oscil- 
lations requires almost maximal mixing. The MSW solu- 
tions include one for small mixing angle (SMA) and one 
with essentially maximal mixing (LMA). 

Another piece of evidence for neutrino oscillations 
is the atmospheric neutrino anomaly, observed by 
Kamiokande, SuperKamiokande, 1MB, MACRO, and 
Soudan-2 (3). Of these, the Superkamiokande data has 
especially high statistics - roughly 52 kton-years worth 
of data at present. This data can be well fit by the infer- 
ence of v^ -^ Yv oscillations with Aot^„„ ~ 3.5 x 10^^ eV^ 
(3) and maximal mixing sin 20„„„ = 1, where v.v = Vx is 
favored. The possibility v^ = ^sterile is disfavored at the 
2a level (5). (The possibility that v^ is predominantly Vg 
is ruled out by both the Superkamiokande data and the 
CHOOZ experiment (4)). 

In addition, the LSND experiment has reported ob- 
serving v^ -^ \e and v^ -^ Vg oscillations with Anij^^j^j^ ^ 



0.1 — 1 eV^ and moderately small mixing angle. This re- 
sult is not confirmed by a similar experiment, KARMEN 
(6). 

There are currently strong efforts to confirm and ex- 
tend the evidence for neutrino oscillations in all of the 
various sectors - solar, atmospheric, and accelerator. 
Some of these are currently running: the Sudbury Neu- 
trino Observatory, SNO, the K2K pioneering long base- 
line experiment between KEK and Kamioka. Others 
are in development and testing phases, such as Borex- 
ino, KamLAND, MINOS, mini-BOONE, and the CERN- 
Gran Sasso program. Among the long baseline neu- 
trino oscillation (LBLNO) experiments, the distances are 
L ~ 250 km for K2K, 730 km for both MINOS, from 
Fermilab to Soudan and the proposed CERN-Gran Sasso 
experiments. The sensitivity of these experiments should 
reach the region Arn^ ^ few x lO^^eV^. Another gen- 
eration of experiments, with even higher sensitivity will 
be required for precision measurements of oscillation pa- 
rameters. One of the physics capabilities of the Next gen- 
eration Nucleon decay and Neutrino detector discussed 
at this NNN99 workshop would be as part of a LBLNO 
experiment. An interesting possibility that is being stud- 
ied intensively is a muon collider or storage ring that 
would serve as a source of quite high intensity, flavor- 
pure (v^, + Ve beams from /j^ and v^ + Vg beam from /j+) 
(anti)neutrino beams. Using these, one could perform 
LBLNO experiments with an existing deep underground 
detector, e.g., at Soudan, Gran Sasso, or Kamioka, the 
NNN detector, and/or a surface detector. Studies have 



shown that one can get hundreds of events per kiloton- 
year at distances of 7000-9000 km (7), (8). It is thus 
appropriate to begin planning for this next generation of 
very long baseline neutrino oscillation experiments. 

An important effect that must be taken into account 
in such experiments is the matter-induced oscillations 
which neutrinos undergo along their flight path through 
the Earth from the source to the detector. In a hypotheti- 
cal world in which there were only two neutrinos, v^ and 
Vx, the v^ -^ Vx oscillations in matter would be the same 
as in vacuum, since both have the same forward scatter- 
ing amplitude, via Z exchange, with matter. However, 
in the realistic case of three generations, because of the 
indirect involvement of Vg due to a nonzero L'^13, and be- 
cause of the fact that Vg has a different forward scattering 
amplitude off of electrons, involving both Z and W ex- 
change, there will be a matter-induced oscillation effect 
on v^ -^ Vx (as well as other channels). An early study of 
matter effects in the earth is (9); several recent studies are 
(8)-(15). 

Here we shall report on a study that we have car- 
ried out (16) of matter effects relevant to LBLNO experi- 
ments. We consider the usual three flavors of active neu- 
trinos, with no light sterile (= electroweak-singlet) neu- 
trinos. This is sufficient to describe the more established 
evidence from the solar and atmospheric neutrino deficit 
(if one tried to fit also the LSND experiment with a neu- 
trino oscillation scenario, one would be led to include 
light sterile neutrinos). As suggested by the solar and 
atmospheric data, we consider that there is only one mass 
scale relevant for long baseline and atmospheric neutrino 
oscillations, Am^,^ ^ few x 10^^ eV^ and we work with 
the hierarchy 
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In our work we take into account the actual profile of 
the Earth, as given by geophysical seismic data (17) and 
compare the results with those calculated using the ap- 
proximation of average density along the path of the neu- 
trino. Further, we present the oscillation probabilities as 
functions of E/Arri^ so one can determine which ener- 
gies are best suited for precise measurements of Am^ in 
a given region. We study how these oscillation probabil- 
ities vary with the different input parameters and discuss 
the influence of the matter effects on the sensitivity to 
each of these parameters. 



MATTER EFFECTS 



The evolution of the flavor eigenstates is given by 
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where the flavor neutrino wavefunction is 

V = UV„r 

in terms of the mass eigenstates 



and 







M^= I mj I ,y = 

mj^ 














V 






(3) 



(4) 



(5) 
where A^^ is the electron number density and 

^/2GFNe[&y] = 7.6 x lO-'-^y^p [glcrc?]. 



The leptonic mixing matrix U can be written as 
U = R2iKRnK*Ri2K' 



(6) 



(2) 



which is the standard CKM-type parametrization, with 
Rij being the rotation matrix in the ij subspace, c\2 = 
cos 012, ^12 ~ sin9i2, etc., K — diag{e^ , 1,1), and K' — 
ii/flg(e'"i,e'"2 1) (the latter phases originate from the 
general presence of Majorana mass terms but will not be 
important here). 

The atmospheric neutrino data suggests almost max- 
imal mixing in the 2 — 3 sector. However, a small but 
non-zero ^13 is still allowed, and this produces the mat- 
ter effect in the traversal of neutrinos through the Earth. 
We use sin^(29i3) < 0.1, consistent with the limits from 
the atmospheric neutrino data (3) and the CHOOZ ex- 
periment (4). We also assume the small mixing angle 
(SMA) MSW solution to the solar neutrino data. This as- 
sumption, together with the hierarchy of eq. (1), implies 
that, for the relevant energies £ > 1 GeV and pathlengths 
L ^ 10^ — 10'* km, only one squared mass scale, Am^(„, is 
important for the oscillations, and the three-species neu- 
trino oscillations can be described in terms of this quan- 
tity, Am^,„j, and the mixing parameters sin (2823), and 
sin (2813), with negligible dependence on sin (2812) and 
5; hence also, CP violation effects would be negligibly 
smafl here, and P{Va -> Vb) = Pi^h -^ V„), P{Va -^ Vh) = 
P{^b -^ Vfl). Although, a priori, CP violation would lead 
to P{Va — ^ Vi) 7^ P{Va ^ Vfc) in vacuum, this inequality 
is true in matter even in the absence of CP violation. 

For our purposes, we recall that the Earth is composed 
of crust, mantle, liquid outer core, and solid inner core, 
together with additional sublayers in the mantle, with 
particularly strong changes in density between the lower 
mantle and outer core. The density profile of the Earth is 
shown in fig. 1 from (17). The core has average density 
Peore = 11.83 glcrcr' and electron fraction Yg^am = 0.466, 
while the mantle has average density Pmantie = 4.66 g/cm^ 
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FIGURE 1. Density profile of the Earth 
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FIGURE 2. 



and Ygmamie — 0.494. If one approximates the density as 
a constant along the neutrino flight path, the evolution 
equation can easily be solved, with well-known results. 
However, when one takes account of the actual variable- 
density situation in the earth, it is necessary to perform 
a numerical integration of the evolution equation, which 
we have done. 



RESULTS AND DISCUSSION 

For long baseline experiments like K2K, Fermilab to 
Minos, and CERN to Gran-Sasso, the neutrino flight path 
only goes through the upper mantle. The density in this 
region is practically constant, and the oscillation prob- 
abilities can easily be calculated. The matter effects 
are small, but possibly detectable for the longer base- 
lines. However, there are several motivations for very 
long baseline experiments, since, with sufficiently high- 
intensity sources, these can be sensitive to quite small val- 
ues of Arn^ and since the matter effects, being larger, can 
amplify certain oscillations and can, in principle, be used 
to get information on the sign of Am^^^. Hence we con- 
centrate here on these very long baseline experiments; for 
these, the neutrino flight path goes through several layers 
of the Earth with different densities, including the lower 
mantle. We show results for L ~ 7330 km, the distance 
from Fermilab to Gran Sasso. We have also performed 
calculations for the Fermilab to SuperKamiokande and 
Fermilab to SLAC path lengths, -- 9120 and 2880 km, re- 
spectively. We calculate the probabilities of oscillation in 
long baseline experiments as a function of E/Arn^, rather 
than using a particular value for Arri^ or the energy. The 
relevant ranges are Am^ ^ few x 10^^ eV^ and energies 



E of the order of tens of GeV. This way of presenting the 
results can be useful in studying the optimization of the 
beam energy. In our work (16) we calcualte the oscil- 
lation probabilities for different values of the mixing an- 
gles 9i3 and 823 allowed by the atmospheric neutrino data 
and the CHOOZ experiment; for this workshop report we 
only show results for sin (2823) ~ 1. We consider both 
neutrinos and antineutrinos. The matter effects change 
sign in these two cases; for antineutrinos, V in (5) is re- 
placed with (—V). This implies that if Am^ is positive 
(as considered here), one can get a resonant enhancement 
of the oscillations for neutrinos, while for antineutrinos 
the matter effects would suppress the oscillations. The 
situation would be reversed if Arri^ were negative. 

We first study the survival probability of v^. If the 
beam went through vacuum, the oscillation probability 
would look like the curve in fig. 2 for practically any 
value of sin (2813). In matter, this probability becomes 
sensitive to all oscillation parameters, as can be seen from 
fig. 3 and fig. 4. 

We also want to compare the solution in vacuum (fig. 
2) with the solution in matter for neutrinos (fig. 3) and 
antineutrinos (fig. 5). In the legends for the figures with 
antineutrinos, "anti Va -^ Vh means Va -^ V/,. One can see 
the opposite effects of matter on neutrinos and antineu- 
trinos. The difference in the results for different mixing 
angles makes it possible in principle to use this proba- 
bility for relatively precise measurements of the oscil- 
lation parameters. Measuring separately the probability 
for V and v can be very useful in detecting the matter ef- 
fects and using these to constrain the relevant mixings and 
squared mass difference. Clearly, if one could use two 
path lengths, as may be possible with a neutrino factory, 
this would provide more information and constraints. 
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FIGURE 5. 



The relative effects of matter can be especially dra- 



matic in the oscillation probability P(v^ 
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If the beam were to go through 
would be probably too small to 



these directly involve v^. 
the vacuum, P(v^ -^ Ve] 
detect (fig. 6). Because of the matter effect however, this 
probability can be strongly enhanced, as is evident in fig. 
7. The enhancement is largest for E/Ain^ around 3000 
GeV/eV^. This is essentially equal to the ratio that one 
would get using a beam energy of ^ 10 GeV, given the 
indication from the data that Aniatm = 3.5 x 10^^ eV^. 
Hence the matter effect can amplify P{Vp — > V^.) and en- 
able this transition to be measured with reasonable accu- 
racy, thereby yielding important information on the os- 
cillation parameters. This probability is very sensitive to 
the value of 9i3 (figs. 7, 8), so one could use it for a 
good determination of this angle. The sensitivity to Am^ 
is also quite strong, due to the pronounced peak given by 
the matter effect in the relevant region. Note that for an- 
tineutrinos, the oscillation is suppressed (fig. 9), so an 
independent measurement of the two channels (v^, ^ Vg 
and v^ -^ Ve) would be very valuable. 

The atmospheric neutrino data tells us that the dom- 
inant oscillation channel is actually v^, — > Vx. Conse- 
quently, it would be very useful to measure P(v^ -^ Vx); 
this would provide further confirmation of this oscillation 
and could also provide accurate determinations of Am^ 
and 623. Fig. 10 shows P(v^ -^ Vx). Fig. 11 shows 

P(v^^Vx). 

Since with a muon collider or muon storage ring, Vg 
(Vj.) beams would also be available, it would be interest- 
ing to study oscillation probabilities with these beams. 
We already have the results for P{Ve ^ v^) since, as men- 
tioned above, with our parameters, this is the same as 
^(v^ -^ Vg). We present here P(Ve ^ Vx) in fig. 12 and 
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FIGURE 10. 
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P{Ve — ^ Vx) in Fig. 13. These calculations show that mat- 
ter effects are important and enhance oscillations of the 
neutrinos and suppress oscillations of antineutrinos in the 
relevant region of parameters. 

To summarize, in planning for very long baseline neu- 
trino oscillation experiments, it is important to take into 
account matter effects. We have performed a careful 
study of these, including realistic density profiles in the 
earth. Matter effects can be useful in amplifying neu- 
trino oscillation signals and helping one to obtain mea- 
surements of mixing parameters and the magnitude and 
sign of Ain~. 
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